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Static VAR Compensators (SVCs) is a Flexible Alternating Current 
Transmission System (FACTS) device that can control the power flow in 
transmission lines by injecting capacitive or inductive current components at 
the midpoint of interconnection line or in load areas. This device is capable 
of minimizing the overall system losses and concurrently improves the 
voltage stability. A line index, namely SVSI becomes indicator for the 
placement of SVC and the parameters of SVCs are tuned by using the multi¬ 
objective evolutionary programming technique, effectively able to control the 
power. The algorithm was tested on IEEE-30 Bus Reliability Test System 
(RTS). Comparative studies were conducted based on the performance of 
SVC in terms of their location and sizing for installations in power system. 
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1. INTRODUCTION 

Nowadays, electrical power system becoming more stressed and difficult to operate due to 
increasing electricity demand and also restrictions on the construction of new lines [1], This causes greater 
losses in the power system network. Besides that, due to constraints, voltage and stability limits work below 
their thermal rating. Furthermore, substitution of new transmission line cost very high to take increased load 
demand on the power system. Hence, to overcome inefficient utilization of transmission system 
infrastructure, power electronics have been developed the flexible AC transmission system (FACTS) devices 

[2], FACTS devices used to control the main parameter, namely, steady state and transient stabilities of a 
complex power system and voltage [3], This parameter influences the AC power transmission [4], Moreover, 
this device act as ideal switch that is controllable and also sustain high levels of tension and power. 
Normally, FACTS device can be divided into two, namely, the Thyristor-Controlled FACTS devices and 
Voltage Sourced Converter (VSC)-based FACTS devices. The Thyristor-Controlled FACTS devices are 
included SVC and TCSC while VSC based FACTS devices are including IPFC and SSSC, shunt devices like 
STATCOM. In power flow studies Thyristor-Controlled FACTS devices act as controllable impedance, but 
VSC-based FACTS devices act as controllable sources which is more complex [5] [6]. 

Recently, researcher has been developed algorithm incorporating with the FACTS device in order to 
solve the problem. The techniques are included. Genetic Algorithm (GA), Bee Algorithm (BA) and Particle 
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Swarm Intelligence (PSO) [7]-[9], In brief, the FACTS device used sensitivity based method to determine the 
optimal placement of FACTS devices for voltage stability improvement. In 2001, Padhy proposed a method 
that incorporates with FACTS to overcome optimal reactive power dispatch problem. The proposed method 
is a new hybrid model for optimal power flow. The fuzzy variables of the proposed method are including 
generation outputs, cost of generation and load demands [6]. Besides that, Chung and Fi suggested an 
improved GA to overcome problems regarding optimal reactive power dispatch by using FACTS device. 
There are two types of FACTS device have been used to control power flow. It is included TCPS and TCSC. 
The developed algorithm that is GA joined with AC power flow in order to reduce the total generation fuel 
cost and to maintain the power flow within the constraints [10]. Similarly, in 2006, Shao and Vittal proposed 
a linear programming (FP) with a FACTS device for solving Optimal reactive power dispatch problem that is 
lightening overloads in transmission line and voltage violations due to contingency. The developed algorithm 
applied to the New England 39-bus system and the WECC 179-bus system [11]. 

This paper presents a new approach for installation of FACTS based SVC on Multi-Objective 
Evolutionary Programming (MOEP) optimization technique considering multi-contingencies (N-m) 
occurrence in the system. The proposed technique determines the optimum sizing of Static VAR 
Compensator (SVC) in order to reduce the total transmission loss in the system. Static Voltage Stability 
Index ( SVSI ) is used as the tool to indicate the SVC’s location to be installed into the power system network. 
The SVSI and transmission loss minimization was used as the objective function in the system. A computer 
program was written in MATFAB and the proposed techniques were tested on the IEEE 30-bus RTS. In 
addition, comparative studies are conducted by comparing the results with Multi-Objective Artificial Immune 
System (MOAIS). An algorithm to apply such multi-objective optimization has been formulated based on the 
same non-dominated sorting concept implemented in non-dominated sorting genetic algorithm (NSGA-II). In 
addition, a program is also developed to obtain best compromise solution in a power system. 

2. MULTI OBJECTIVE OPTIMIZATION 

Multi-objective optimization is a process to find the value of the variables that minimize the 
objective function namely SVSI and transmission loss while the system is operating within its constraint limit. 
Multi-objective problems are more difficult to solve compared to the single objective since there is no unique 
solution. Instead of one optimal resolution, the implementation of multi-objective can give a set of optimal 
solutions. These optimal solutions are known as Pareto-optimal solutions. The set of all feasible non- 
dominated solution is referred to as the Pareto optimal set, and for a given Pareto optimal set, the 
corresponding objective function values in the objective space is called the Pareto front. The multi objective 
optimization problem is specified as follows [12]: 

Min / maximization : F(x) = [fi(x), f 2 (x), f 3 (x) . .f k (x)] 

Subjected to : g, (x)< 0 where i= 1,2,3 . i (1) 

: hj(x) = 0 where k= 1,2,3. .j 

Where F(x) is objective, ffx), f 3 (x), f 3 (x) . f k (x) are the objective function, x is the vector of 

dependent variable, g is the equality constraints and h is the inequality constraint. 

The FACTS device installed on the weak buses and heavily loaded areas in order to reduce the 
stress condition in the system. The locations of SVCs devices indicated using the SVSI technique that 
operates at same operating conditions in the power system network. When the load flow program was run, 
stability indices are calculated and the system identified the lin with the highest SVSI fo the installation of 
FACTS device. Finally, MOEP technique was used to identify the optimal size of the SVC. The process of 
installation considered the occurrence of generator outages. 

2.1 Static Voltage Stability Index ( SVSI ) 

SVSI which is a line-based voltage stability index was developed by [13]. This index used in the 
voltage stability analysis as an indicator of the voltage stability condition of a system. The voltage stability 
condition of all lines in power system could be assessed using this index, which could predict the occurrence 
of voltage collapse in a system. SVSI was formulated by deriving the voltage quadratic equation for a general 
two-bus system at the receiving end. SVSf, can be defined as shown in Equation (2) for the two-bus system. 
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Where the active power and reactive power are Pj, and Q jr The line resistance and reactance are R Jt 
and Xj h The voltage magnitude and angle are IVI and S. The subscript i and j denote variables associated with 
bus i and bus j. SVSI indicated the steady state voltage stability of the line. If the SVSI is less than one, there 
are solutions and the system is stable. If the SVSI is larger than one, there is no solution and the system 
becomes unstable or steady state voltage collapse occurs in the system. 

2.2 Minimization of transmission loss as objective function 

Another objective function considered in the proposed method is minimizing the transmission power 
losses in the transmission network, while satisfying a set of physical and operation, subjected to a set of 
equality and inequality constraints in the power system [14], 

The mathematical equation of transmission loss can be written as 

min 4= IX„M)= ^g k {yr+Vi-wyj cos 0^) 

k<sN E ksN E 

M'-.d (3) 

Subject to:- 

Ki = Qa - Qd, - V t Y, V ji G :j sine n ~ B u cos = n pq 

je N, 

Subject to the constraint of equality in reactive and active power balance 

Q, ~ Qa + Qm = 0 . Q, = Qa ~ Qa ~ V i T/j i G iJ sin d U ~ B U cos B «) = 0 ’ ' e N r Q 

m, (4) 

P-P Gi +P Di = 0, P i =P ci -P Di -V i Yy j (G ij cose ij -B ij vnd i )=0, isN^ 

jcN, 

Hence, inequality constraints on control variable limits; generator power reactive capability limits, 
generator power active capability limits, and voltage constraints are given by; 

2c, ni „ ^ Qa < Qa . > e N g 

Q cl .^ Q ci ^ Q d _ i e V ( 5 ) 

P Gi <P Gi < P Gj , i e Slackbus 
V ijnh <V,< V L _ i e N b 

where, g k is the conductance of branch k, n s is the slack (reference) bus number; N P q is PQ bus 
number, Npv is PV bus number, N B is the total number of buses, N P -i is the total buses excluding slack bus, N c 
is the possible reactive power source installation buses number, N l: is the branch number, N, is the numbers of 
buses adjacent to bus i including bus i, 0 :j is voltage angle different between bus i and bus /(rad), Q, and Q, 
are the reactive power on the sending and receiving buses; Q c is the generated reactive power, V, and V, are 
the voltage magnitude at the sending and receiving buses , G,j and B tJ is the mutual conductance and 
subceptance between bus i and bus j and P K is the total active power loss in the system. 


3. MULTI-OBJECTIVE EVOLUTIONARY PROGRAMMING 

The MOEP mainly carried out six steps, namely, initialization, non-dominated sorting, crowding 
distance, mutation, combination and selection. The population is initialized to generate random number 
generation. In MOEP initialization is one of the important processes to produce first population termed as 
parents. Then, The population is sorted based on the non-domination. Each solution should be compared with 
every other solution in the population to find if it is dominated. Each solution assigned a fitness or rank equal 
to its non-domination level (1 is the best level, 2 is the next best level and so on). Furthermore, the first rank 
belongs to the most excellent non-dominated set in the population [15]. Once the non-domination sort is 
completed, the crowding distance is assigned. Crowding Distance also known as a fitness value of an 
individual. The purpose of crowding distance is to provide the diversity in the population [12]. Then, the 
individual solutions are sorted in descending order based on the magnitude of the crowding distance values. 
Subsequently, the process continued with the mutation process. The mutation operator changed its current 
value of a continuous variable to a neighboring value using Polynomial Probability Distribution and this is a 
basic procedure of any genetic operator [16]. The offspring produces from the mutation process are combined 
with the clone parent to undergo a selection process in order to identify the candidates have the chance to be 
transcribed in the following generation. The best individual from the offspring population will be selected 
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according to a selection scheme in order to form the parent population for the following generation. The 
offspring is chosen using tournament scheme. Finally, the MOEP produces a set of Pareto optimal solution 
which one objective cannot be improved without sacrificing other objective. Therefore, from the Pareto- 
optimal set of non-dominated solutions, the Best Compromise Solution (BCS) was selected for the decision 
maker as the BCS will decide the best solution in between both objective functions [17]. 



O.-v-Qmm ~ Q.\ht 


Figure 1. SVC Model 


4. STATIC VAR COMPENSATOR (SVC) 

Nowadays, SVC is the most popular FACTS devices which are used to solve optimal reactive power 
problem (ORPD). Furthermore, SVC not only generate reactive power but is also absorbing reactive power. 
SVC connected in parallel to transmission line where TCR located in parallel with capacitor bank [18]. 
Usually, SVC installed at the end of the transmission line or midpoint of transmission interconnections. 
Moreover, SVC is a three phase and shunt connected device. The main function of SVC is to improved 
voltage in weak transmission line. 

The mathematical modelling of SVC considered in this study is shown in Figure 1. In recent years 
SVC has been used for reactive power support and to enhance voltage stability in the electrical power system 
network. Furthermore, the SVC can act as both inductive and capacitive compensation by absorbing releasing 
reactive power. Hence, it is modelled as ideal reactive power injections to perform the steady-state condition 
at bus i. The absorbed or injected power at bus i in the system is represented by Q svc . The mathematical 
formulation of SVC constraint shown as 


Qmin — Qsvc - Q max 

-200 MVar < Q svc < 200 MVar 


5. APPLICATION OF MOEP IN SVC DEVICE INSTALLATION 

MOEP involved initialization, non-dominated sorting, crowding distance, mutation, combination 
and selection. Transmission loss minimization and voltage stability were chosen as the objective function for 
the optimization process. The flow chart for the application of MOEP in SVC device installation is shown in 
Error! Reference source not found.. Several inequality constraints are set in this study so as to achieve the 
optimal solution. Several generator outages namely generator at bus 11 and 13 were subjected into the 
system. The selections of outages are based on the most severe generator and in the system to maximize the 
performance of the system. There are two constraints assigned before the SVCs sizing is optimised. The 
constraints are; total loss to be less than the loss_set and voltage at the loaded bus higher than V_set. The 
loss_set and V_set are the total loss and voltage at the loaded bus before the optimisation process was 
conducted. The following steps show the implementation of EP. 

i. Set the generator outages. 

ii. Set the loading factor, 1. 

iii. Setting the location for SVC using SVSI index of stability. 

iv. Set the ORPD constraints i.e. SVSI < SVSI_set and total loss < loss_set as objective functions. 

v. Generate random number i.e. x h x 2 ,... x 5 . Check for constraint violations. If constraints violated, 
go to step iv, otherwise go to step vi. 

vi. Fill in population pool. Repeat step (ii) if pool was not full, otherwise continue to step (vii). 

vii. Calculate the non-dominated solution for each individual in the population. 

viii. Sort the entire population using front. 

ix. Calculate the crowding distance for each front 
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x. Mutate the parents x h x 2 , x 3 , x 4 and x s to generate offsprings. 

xi. Assign x h x 2 , x 3 , x 4 and x 5 in the bus system data. 



Set the generator outages (N- 
m) 


Set the loading factor 


Run load flow 


Calculate index 


Sort index in descending order 


Displays the highest index 
(Sensitive Line) 


Set the constraints 


Generate random no, x as a 
control variable (x],x 2 ,...x n ) 
n= no. of installation 



Calculate the non-dominated 
front for each individual in 
population 


Sort the entire population 
using front 


Calculate the crowding 
distance for each front 


Mutate the parents 
(Generate offsprings) 


Assign x to system 


Calculate fitness by running 
load flow programme to 
evaluate SVSI and total loss 


Combine parent and 
offsprings 

I 

Perform selection by ranking 
process 




Yes 

Sort the Pareto optimal front 



Find best compromise solution 



Plot the Pareto optimal front 
=1 and BCS into graph 




Figure 2. Flowchart for SVC installation using MOEP 


xii. Calculate fitness by running load flow program to evaluate SVSI values and transmission loss 
values. 

xiii. Combine parents and offspring (combination process). 

xiv. Perform selection by tournament selection process from the combine data. 

xv. Identify and transcribe new generations. 

xvi. If solution is not converged, repeat step v to xii, otherwise go to step xvii. 

xvii. Sort the Pareto optimal front, 
xviii. Find the best compromise solution. 

xix. Plot the Pareto optimal front and the best compromise solution into graph. 

xx. Stop 


6. RESULTS AND ANALYSIS 

The results have been obtained from the developed algorithm for multi-objective evolutionary 
programming for SVC installation and sizing. The developed algorithm has been tested on the IEEE 30-Bus 
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RTS. The SVC installation has been formed by considering multi-contingencies ( N-m ); generator outages 
have been occurred in the system which are formulated with the objective of minimizing the SVSI and real 
power loss, respectively. For the studies, the population size selected is 200, mutation probability is 0.1 and 
distribution index for mutation is 20. Population size of 200 was chosen based on exhaustive experiment on 
various population sizes which revealed that 200 is the most suitable one. In this study, MOEP algorithm has 
been applied where both SVSI and transmission loss were optimized simultaneously. The developed 
algorithm has been tested with bus 26 subjected to a maximum loadability in the IEEE 30-bus RTS. The 
simulation results are tabulated in Table 1. The Pareto optimal fronts set are shown in Figure 2. 

The diversity of the Pareto-optimal set over the trade-off surface for bus 26 is shown in this figure. 
It is observed that the proposed approach produces 200 Pareto-optimal solutions. Out of the solutions 
obtained, two non-dominated solutions that represent the best SVSI value and best transmission loss are given 
in Table 1. The best compromise solution of this method is also given in this table. It is observed that the 
SVSI and transmission loss values, reduce with respect to 1 after the implementation of MOEP in the system. 



Figure 3. Pareto front for SVSI and transmission loss minimization obtained using MOEP for SVC sizing at 

bus 26 


Table 1. Results for MOEP when bus 26 wa s reactivelv loaded for \= 30 MVAr 


Units 
(in MY Ar) 

pre-SVC 

Solution at niinininni 

srs/ 

MOEP 

T ransmission Loss 

Best Compromise 
Solution 

SVC, 

Bus 26 

-3.2409 

-1.8782 

-1.9755 

SVC; 

Bus 27 

-32.2077 

-30.0446 

-31.0392 

SVC, 

Bus 25 

20.0739 

19.5343 

19.0337 

SVC, 

Bus 27 

-32.5618 

-33.4772 

-33.2473 

SVC, 

Bus 5 

-27.4607 

-25.4652 

-27.8897 

.sits/ 

0.5835 

0.16787 

0.16792 

0.167S9 

Transmission Loss 
(in MW) 

28.0137 

16.8913 

16.8777 

16.8794 


Table 2. 

Results forMOAIS whenbus26 was re actively loaded for /.= 

30 MYAr 

Units 
(in MYAr) 

pre-SYC 

■Solution at minimum 

srs/ 

MO.VIS 

Solution at min ini nm 

T ranszuission Loss 

Best Compromise 
Solution 

SVC, 

Bus 26 

-9.0485 

-6.6835 

-7.0939 

SVC, 

Bus 27 

-25.7690 

-25.4146 

-25.7839 

SVC, 

Bus 25 

10.7274 

12.2215 

11.9964 

SVC, 

Bus 27 

28.5779 

27.3093 

27.2783 

SVC, 

Bus 5 

-26.5534 

-27.0564 

-27.1382 

SI’S/ 

0.5835 

0.16790 

0.16792 

0.16791 

Transmission Los: 
(in '.TV 

s 28.0137 

16.9335 

16.9031 

16.9066 


Table 3. Best Com 

promise Solution for SVC installation when bus 26 was reactivelv loaded. 

T echnique 

Non dominated 
Solutions No. 

SIS/ 

T ranszuission Loss 
(in MAY) 

MOEP 

196 

0.16789 

16.8794 

MOAIS 

12 

0.16791 

16.9066 
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It implies that the voltage stability has been improved. The best SVSI value is 0.16787 while the best 
transmission loss is 16.8777 MW. The best compromise solution is 0.16789 for SVSI and 16.8794 MW for 
transmission loss. It is observed that the SVSI and transmission loss have been reduced with the 
implementation of MOAIS. It implies that the voltage stability has been improved. The best SVSI value is 
0.16790 while the best transmission loss is 16.9031 MW. The best compromise solution is 0.16791 for SVSI 
and 16.9066 MW for transmission loss. 

The comparison result for the best compromise solution for different optimization technique using 
MOEP and MOAIS for the implementation of SVC is tabulated in Table 3. In this table, the results are 
verified from three aspects in terms of SVSI value, transmission loss and amount of non-dominated solutions. 
When the load is subjected to bus 26, it shows that only 12 non dominated solutions distributed along Pareto 
Front using MOAIS. Nevertheless, the MOEP has presented 196 non-dominated solutions along the Pareto 
Front, in which it gives more choices of selection for SVC installation size to improve SVSI and reduce the 
transmission loss in the system. As highlighted in the table, it is observed that MOEP outperformed MOAIS 
since MOEP managed to improve the SVSI value and transmission losses as compared to MOAIS in the 
system. 


7. CONCLUSION 

This chapter has presented multi-objective optimization techniques termed as MOEP in 
implementing the optimal sizing of the SVC installation scheme. The combination of SVSI and transmission 
loss minimization as objective function has been solved for the IEEE 30-bus RTS system with bus 26 
subjected to loading condition. The Pareto-optimal front has been obtained in all schemes and the best 
compromise solution shows the promising results where MOEP and MOAIS successfully improved the SVSI 
value and reduced the transmission loss values in the system. Finally, the results obtained from the MOEP 
techniques were compared with MOAIS and it was found that MOEP outperformed MOAIS in all cases. 
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